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ABSTRACT 

We present the first results of an ongoing campaign using the STIS spectro- 
graph on-board the Hubble Space Telescope (HST) whose primary goal is the 
study of near ultraviolet (UV) spectra of local Type la supernovae (SNe la). Us- 
ing events identified by the Palomar Transient Factory and subsequently verified 
by ground-based spectroscopy, we demonstrate the ability to locate and classify 
SNe la as early as 16 days prior to maximum light. This enables us to trigger 
HST in a non- disruptive mode to obtain near UV spectra within a few days of 
maximum light for comparison with earlier equivalent ground-based spectroscopic 
campaigns conducted at intermediate redshifts, z ~ 0.5. We analyze the spectra 
of 12 Type la supernovae located in the Hubble flow with 0.01 < z < 0.08. Al- 
though a fraction of our eventual sample, these data, together with archival data, 
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already provide a substantial advance over that previously available. Restricting 
samples to those of similar phase and stretch, the mean UV spectrum agrees 
reasonably closely with that at intermediate redshift, although some differences 
are found in the metallic absorption features. A larger sample will determine 
whether these differences reflect possible biases or are a genuine evolutionary ef- 
fect. Significantly, the wavelength-dependent dispersion, which is larger in the 
UV, follows similar trends to that observed at intermediate redshift and is driven, 
in part, by differences in the various metallic features. While the origin of the 
UV dispersion remains uncertain, our comparison suggests that it may reflect 
compositional variations amongst our sample rather than being predominantly 
an evolutionary effect. 

Subject headings: cosmological parameters — supernovae: general — ultraviolet: 
general 



1. INTRODUCTION 



Type la supernovae (SNe la) remain the most practical and well-exploited cosmological 
probe offering an immediate route to understanding 'dark energy'. Measures of distant 
events are being used to distinguish between Einstei n's cosmological constant, A, and a 
scalar field whose equation of state parame ter w 7^— 1 (lAstier et al.ll2006t iRiess et al.l 120071 ; 
Kessler et al.ll2009l ; lAmanullah et al. II2010I ). Yet despite remarkable observational progress, 
there is no satisfactory theory explaining a SN la event. The mechanism by which a white 
dwar f accretes additional material is unclear as is the nature of the explosion itself (ILivio 
1999h . 



To facilitate progress, observers employ a variety of empirical correlations to reduce the 
intrinsic scatter of the SN la Hubble diagram. SNe la were initially considered a one or 
two parameter family with the light curve width and rest-frame color as the key variables. 
However improved data have revealed important correlations with the host galaxy. Events 
are not only more common in star-forming hosts per u nit stellar mass but their light curve 
properties differ from those seen in quiescent galaxies ( Sullivan et al. 2006. 12010 ). an effect 



that has direct consequences for their use over large look-back times (IHowell et al.l 120071 ). 



These discoveries naturally raise the question of what further evolutionary changes might 
be present in the SN la population. A long-standing concern has been the unknown effect of 
an evolving progenitor composition, both in terms of a possible redshift-dependent bias, and 
in producing an intrinsic dispersion that could limit the effectiveness of large future surveys. 
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A Keck study of high quality rest-frame near-UV SNe la spectra at intermediate redshift 
(z ~ 0.5) reveals a surpr ising diversity a t short wavelengths where some models predict a 
sensitivity to metallicity (jEllis et al.ll2008l hereafter E08, see also Foley et al. 2008 ). A large 



U-ban d dispersion had earlier been claimed in the local photometric survey of iJha et al. 



( 2006 ). although subsequent photometry has challenged the amount ( Astier et al. 20061 ). 
Although so me models predict a UV dispersion mig ht arise from variations in the progenitor 
composition (ILentz et al.ll2000t iHoflich et al.ll2000l ). the magnitud e of the effect seen b y E08 
exceeds that predicted for reasonable compositional differences. ISauer et al. I ( 120081 ) have 
shown a large fraction of the UV flux can be formed by reverse-fluorescence scattering which 
affects the dependence on composition. If the observed dispersion is found to arise from some 
evolutionary trend, it c ould bias futur e z > 1 SNe campaigns that typically sample from 
this wavelength region. ISullivan et al.l (120091 hereafter S09) compare mean SN la spectra 
over a redshift path of < z < 1.2 and find no strong evolution, however, only three local 
UV spectra were available at the time, seriously limiting the comparison and giving no local 
measure of the UV dispersion. 

Efforts to understand the UV behavior of SNe la have subse quen tly intensifie d . Us- 
ing the UV optical telescope on-board SWIFT, brown et all fcoioh and iMilne et al.l (boioh 
have confirmed the presence of a dispersion increase to shorter wavelengths in local SNe 
la. Comparisons between local and intermediate redshift data remain confused and so the 
question remains as to whether some component of the significant diversity of UV spectra 
seen in distant SNe la is an evolutionary phenomenon or represents some as yet unexplained 
diversity in the SN la mechanism. 

As described in E08, spectroscopic studies offer a major advantage over photometric 
investigations as they eliminate uncertainties arising from k-corrections and, with adequate 
data, the mean and dispersion can be investigated in the context of known metallic features. 
Following the successful repair of the UV-capable Space Telescope Imaging Spectrograph 
(STIS) aboard the Hubble Space Telescope (HST) during the 2009 Servicing Mission 4, it 
has become possible to make significant progress in addressing the above questions. Using 
a non-disruptive Target of Opportunity campaign (GO 11721, PI: Ellis), we are securing 
maximum light STIS spectra for 35 SNe la located in the Hubble flow. Here we present the 
results from the first 12 events from this program, augmented by three events from earlier 
archival data, that provide comparable statistics to the survey of E08 at intermediate redshift 
for the analysis undertaken here. 
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Table 1. Low Redshift Supernova la Sample 



Supernova 


Phase a 


Stretch 


Host z 


PTF-09dlc 


+2.02 ±0.38 


1.152 + 0.054 


0.0675 


PTF-09dnl 


+0.68 ± 1.99 




0.0231 


PTF-09dnp 


+5.47 + 0.87 


0.996 + 0.308 


0.0373 c ' d 


PTF-09fox 


+ 1.66 + 0.40 


1.016 + 0.109 


0.0718 


PTF-09foz 


+2.59 + 0.38 


0.883 + 0.079 


0.0543 c ' d 


PTF-lObjs 


+2.10 + 0.18 


1.138 + 0.019 


0.0300 c ' d ' e 


PTF-lOfps 


+8.15 + 4.85 


0.980 + 0.530 


0.0215 c ' d 


PTF-lOhdv 


+3.10 + 0.36 


1.077 + 0.064 


0.0533 


PTF-lOhmv 


+2.73 + 0.09 


1.150 + 0.009 


0.0324 


PTF-lOicb 


+0.76 + 0.13 


1.071 +0.021 


0.0086 c ' d 


PTF-lOmwb 


-0.19 + 0.14 


0.896 + 0.018 


0.0313 


SN 20091e 


-0.32 + 1.99 




0.0178 d 


SN 1981b 


+1.30 + 0.14 


0.89 + 0.02 


0.0060 f 


SN 1992a 


+5.45 + 0.04 


0.82 + 0.01 


0.0061 f 


SN 2001ba 


+4.21+0.22 


1.02 + 0.02 


0.0305 f 



a Effective phase (phase/stretch) at the time of the 
STIS observations 

b Unless otherwise noted, redshifts are obtained from 
the host features in the PTF spectra and are accurate 
to z < 0.001 

c Sloan Digital Sky Survey (SDSS) 

d NASA Extragalactic Database (NED) 



e May reside in a possible line-of-sight dwarf galaxy 
at z = 0.019 



f Archival UV speci 


ra discussed by S09 ( 


Branch et al. 


1983; 


Kirshner et al. 


1993; 


Folev et al. 


2008) 
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OBSERVATIONS 



A significant challenge in delivering targets to HST for observations at maximum light is 
the need to detect and identify convincing SNe la candidates soon after explosion. For a non- 
disruptiveTarget of Opportunity (ToO) campaign, defined as one where Phase II observations 
are submitted for inclusion in the HST schedule built for the second week following the 
submission, SN confirmation must occur ~ 10 — 12 calendar days prior to submission. As a 
result, spectroscopic SNe la target confirmations ~ 7 — 16 days before maximum light are 
necessary to acquire STIS observations with phases ±4 days from maximum light and match a 
distribution of intermediate- redshift near UV spectra secured by the E08 campaign. Figured] 
illustrates the timeline of the data acquisition for the 12 SNe of our HST program presented 
here - from detection, to ground-based spectroscopic confirmation, to STIS observation - 
and shows the ability of PTF to detect SN la outbursts as early as 16 days prior to maximum 
light. 

The identification of SNe la suitable for HST non- disruptive observations proceeded in 
two stages. The early identification of candidate events was based on photometric survey 
data. Eleven of the 12 newly-discovered SNe were ident ified between 2009 August 19 an d 
2010 June 29 by the Palomar Transient Factory (P TF; 



Rau et al 



2009J; lLaw et all 120091 ). 



As with the CFHT Supernova Legacy Survey (SNLS lAstier et al.ll2006l ) utilized by the E08 
intermediate redshift program, PTF is a rolling search for transient events unbiased to the 
nature of the host galaxy. The re maining event (SN 2 0091e) was triggered on the reported 
results of an independent search (jPignata et al.l 12009k CHASE) during a period when the 
Palomar observatory was under extensive ash clouds from the 2009 California forest fires. 

Ground-based spectroscopic follow-up acquired within one to a few days after photo- 
metric detection represents the second stage of confirmation essential for determining the 
SN type and redshift. Spectroscopy was performed using regularly scheduled time and occa- 
sional ToO interrupt s on the following telescopes/instrume nts: the Low-Resolution Imaging 



Spectrometer (LRIS; l0ke et al. 



Spectrograph (DEIMOS; 
Spect rograph (GMOS; 



19981 ) and X-shooter (IVernet et al 



Spectrograph (DBSP; lOke fc Gunn 



1995 



Faber et al. 



Hook et al 



McCarthy et al.lll998l) and DEep Imaging Multi-Object 



20031 ) on the Keck telescopes, the Ge mini Multi-Object 



2004h on the Gemini telescopes, FORS (lAppenzeller et al.l 



2009) on the ESO Very Large Telescopes, the Double 



19821 ) on the Palomar Hale Telescope, and the Interme- 



diate dispersion Spectrograph and Imaging System (ISIS) on the William Herschel Telescope. 
By matching flux- calibrated spectra with selections from an extensive spectral database 
(IHowell et al.ll2005l ). improved phases typically accurate to ±2 days were obtained. 



Near-UV spectra were acquired using the Space Telescope Imaging Spectrograph (STIS; 
Woodgate etaU Il998l : iKimble et all Il998h between 2009 September 02 and 2010 July 13 
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Fig. 1. — Timeline for the HST non-disruptive ToO program (GO 11721, PI: Ellis). The magnitude 
and phase (number of days relative to maximum light) are shown for the PTF, and one non-PTF, 
photometric discoveries (squares), ground-based spectroscopic confirmations (triangles), and near- 
UV STIS spectra (diamonds) . Photometric discoveries and spectroscopic confirmation are necessary 
during the time windows indicated by the light-gray and gray block regions, respectively, in order 
for STIS near UV spectra to be acquired within the time window indicated by the dark-gray region. 
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Fig. 2. — STIS near UV maximum-light spectra as itemized in Table [T] reduced to rest-frame 
wavelengths. Each label gives the PTF identification (with 091e indicating SN 20091e) and phase of 
STIS observations. Spectra in red represent those not included in the mean and dispersion analysis 
(see text). 
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(Fig. [U red diamonds). To best match our z ~ 0.5 Keck spectra (E08), we used the 
STIS/CCD 430L observing mode with 3 — 4 CR-SPLIT exposures totaling 1 orbit, ensuring 
good signal-to-noise ratio coverage down to a rest wavelength below 2900A. A montage of 
the STIS spectra is shown in Figure [2J Associated near-simultaneous ground-based optical 
spectra were also taken using further regularly scheduled time and ToO interrupts using 
the facilities described above. All data were reduced using standard IRAF and IDL data 
reduction routines. Pre- and post-maximum light PTF r-band photometric observations 
were used to determine th e initial phase and s tretch for each HST triggered event using the 
SiFTO light curve fitter (jConley et al.ll2008l ). Details for the SN la sample here and the 
three archival SNe la (discussed in S09) are listed in Table [TJ 

In 4 cases (09dnl, 09dnp, 09fox, lOfps), the PTF light curves are of marginal quality or 
too poorly sampled for an accurate determination of one or both of the phase and stretch. 
One event (SN 20 091e) has no PTF phot ometry and its phase was determined from the 
CHASE estimate ( iChallis fc Berlindl 120091 ) which is consistent with our spectroscopic esti- 
mate. In our final analysis (§3), this imprecision affects the discussion of two of our STIS 
spectra (09dnl, 20091e). 



3. ANALYSIS 

Prior to HST Cycle 13, UV spectra were secured for 13 local SNe la, but only 5 were 
studied close to maximum light and 2 of these were peculiar. Those suitable from earlier 
work are listed in Table [1] and shown in Figure [2j In Cycle 13 a ToO campaign began 
to increase the sample (GO 10182, PI: Filippenko) but the failure of STIS curtailed this 
program. Remarkably, prior to our study, more was known about the pre-maximum and 
maximum light UV spectra of SNe la at z ~ 0.5 than at z — 0. Our STIS sample has now 
dramatically improved this situation. 

In order to compare the mean UV spectra of local SNe la and their dispersion with the 
sample discussed by E08, we must match the phase (and ideally the stretch) distributions 
of the two samples. Such a matched comparison was not possible in S09 due to the paucity 
of the local data. To maximize the utility of the new HST data, we adopt a phase range of 
—0.32 to +4d (Figured]) bounded by the earliest phase of the HST data (—0.32) and a +4d 
phase limit similar to that applied in E08 to minimize phase evolution effects on the mean 
and dispersion. This criterion leads to 11 events from the HST+archival sample drawn from 
Tableland 16 events from E08. Ignoring the two events for which accurate stretches cannot 
be determined, the mean values are Shst = 1-030 ± 0.038 and S z ~o.s = 1-049 ± 0.019 which 
are consistent within the errors. Excluding the two supernovae for which reliable stretches 
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Fig. 3. — Mean near-UV spectra for low- and intermediate-redshift SN la near maximum light. 
The mean near-UV spectrum derived from the 10 STIS supernovae and one archival SN compliant 
with the adopted phase and stretch criteria (see text) is shown by the blue line with the region 
containing 90% of the jack-knife resampling fits shown as the light blue region. This is compared 
to the mean of a z — 0.5 E08 sample of 16 SNe closely matched in phase and stretch (red line) and 
its 90% region (bounded in pink). 
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are not available does not significantly affect our results in the following subsections. 



3.1. Mean Type la Spectrum 

S09 presented a comparison of t he mean UV spect ra determined by E08 at z ~ 0.5, 



those secured using the ACS grism by lRiess et al.l (12004 ) at z ~ 1.2, and three local spectra 



from archival data listed in Table [TJ Examining the z ~ 0.5 UV spectra, S09 found some 
decrease with redshift in the strength of intermediate mass element features (Si II, Ca II, 
and Mg il), but it was argued this could arise in part due to the natural drift to luminous, 
larger stretch events expected at high redshift. Below A ~ 3600 A the mean local spectrum 
was highly uncertain. 

To facilitate a proper comparison to S09, we construct our mean near-UV spectrum 
following the procedure discussed in E08. Briefly, the spectra are normalized to have the same 
flux through a box filter defined between rest-frame 4000 — 4800A and the variation in the 
mean spectrum is estimated via bootstrap-resampling. Use of other box filters, including the 
full wavelength range common to all spectra, does not significantly affect the results. Figure [3] 
shows the normalized mean spectrum and region containing 90% of 100 bootstrap-resampled 
mean spectra for both the low- and intermediate-redshift matched samples. Although the low 
redshift SN la mean spectrum closely resembles that at z ~ 0.5, as in S09, we notice a marked 
decrease with increasing redshift in the depth of the Si II and Ca II blend near 3800 A, Mg II 
near 4300 A, and those of iron group elements below 3500A. The utility of the comparison is 
more advantageous than that conducted by S09 for wavelengths below 3500A since our new 
local sample is much larger than that used by S09 and the phase and stretch distributions 
are better matched. Farther into the UV we note that the differences between the mean local 
and z ~0.5 spectra become particularly significant. Since the stretch distributions of the 
two samples are similar, this could be a gen uine effect rather tha n one arising from samples 



biased to more luminous and bluer events (IBalland et al.l 120091 ). Analysis of our eventual 
full sample will clarify this important point. 



3.2. Spectral Dispersion 

E08 demonstrated a significant increase in the variance of their z ~0.5 SNe la UV spectra 
for wavelengths below 3700A, both by comparing individual deviations of 15 maximum light 
spectra from their mean in units of the dispersion a, and via photometric colors measured 
directly from color-corrected Keck spectra. For the first time, our low-redshift sample of 
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Fig. 4. — Dispersion from the mean for the 11 low-redshift SNe meeting our phase criteria (multi- 
colored thin curves). The absolute value of the mean dispersion is indicated by the thick black curve 
with la observational uncertainties overlaid (gray region). For direct comparison, the absolute value 
of the mean dispersion for the matched sample of 16 intermediate redshift E08 SNe is shown (thick 
red curve). An increased dispersion at shorter wavelengths is present in both samples. 



-12- 



UV SN la spectra is large enough to perform a similar analysis. Ideally one would color- 
corre ct the local spec tra following the procedure discussed in E08 using the SALT2 color 



law (IGuy et al.l 120071 ) but this requires host- corrected multi-color data that must await 
late-time reference images in bands other than the PTF r-band. We experimented with 
estimating the host contamination in g and i from the contemporary data but concluded the 
uncertainties are too great at this time. Accordingly, a comparison with data from E08 that 
is not color-corrected is more appropriate. It is important to note that the color correction 
E08 applied did not significantly change the UV scatter and its wavelength dependent trend. 
The key question we seek to address is whether the dispersion trend is generic to all SNe la, 
independent of redshift, or largely a feature of the intermediate redshift data only, possibly 
implying some evolutionary effect. 

Figure @] shows that the wavelength-depen dent dispersion is in deed present in the loca l 



data, consistent with the photometric claims of iBrown et al.l ( 120101 ) and lMilne et al.l ( 120101 ). 
In the region that contains 90% of the mean values from bootstrap resampling, the variation 
from one spectrum to another with respect to the mean increases below 3700A as in E08. 
Clearly several of the features which vary between the local and E08 samples discussed earlier 
contribute to the dispersion suggesting a compositional origin. 

Even if we exclude the region dominated by strong features and consider the average 
dispersion from the mean spectrum in the regions UV: 2900-3500 A and optical: 4100-5200 A 
allowing for the spectrophotometric uncertainties, we find that the dispersion increases from 
the optical to the UV by a similar factor of ~2.3 in both the local and z ~0.5 samples. There 
is marginal evidence that the UV scatter may be larger in the z ~0.5 data than in the local 
sample, but the spectrophotometric uncertainties for the individual spectra are naturally 
larger. Overall, we conclude that the UV spectral dispersion is most likely a feature of the 
SNe la population and not an evolutionary effect. 



4. DISCUSSION 

Our initial results clarify and quantify indications from earlier work. Although the SN 
la mean spectrum close to maximum light appears to have remained remarkably similar over 
the past 5 Gyr, we find the decrease in the strength of the metallic features with increasing 
redshift noted by S09 also present in our more representative comparison. Given that the 
mean stretches of the local and z~0.5 samples are similar, this may represent the expected 
decrease in metallicity over this epoch. We can address this possibility in more detail with 
the completion of our survey. 
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Equally important is that we observe a strong wavelength dependent scatter in the 
rest-frame UV spectra of our local sample, as noted in E08. Independent of the calibration 
questions that have plagued recent photometric studies, the spectra demonstrate that the UV 
scatter is generic to SNe la over a wide range of cosmic time and is not likely an evolutionary 
effect. Much of this behavior can be attributed to the varying absorption line strengths of 
intermediate mass elements occupying the UV wavelength region, supporting the notion that 
the UV scatter arises from compositional differences between events. 

In addition to strengthening these conclusions with a larger sample, it is now clear that 
further progress will follow more detailed multi-phase UV studies of selected local events. 
HST has recently been awarded for such a program (GO 12298, PI: Ellis). 

RSE acknowledges support from DOE grant Ide-scOOOllOll MS from the Royal Society, 
A.G. from the Israeli Science Foundation and a European Union Marie Curie fellowship. 
Support for program GO 11721 was provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by AURA Inc. under NASA Contract NAS5- 
26555. 
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